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OPE constraints are studied as a means of distinguishing between the versions of the 1 = 1 vector spectral 
function extracted from (i) inclusive 1=1 hadronic electroproduction cross-sections and (ii) inclusive 7=1 
hadronic t decay data, with the goal of clarifying expectations for the leading order hadronic contribution to 
d/j. = {g — 2);j/2 in the Standard Model. The constraints are shown to, at present, favor the r decay data, 
and hence a Standard Model prediction for compatible with the BNL E821 experimental result. The relative 
role of the tttv and An contributions to the discrepancy between the integrated electroproduction results and 
the corresponding OPE constraints is also investigated, and the significance, in this context, of forthcoming 
re-measurements of the e'^e~ -k^tt'-k^tt^ cross-sections pointed out. 



1. Background 

It is well known that 
• the leading order (LO) hadronic contribu- 



ihad^LO 



tion to a^, [a^ 
dispersive integral 

2 

r ihad.LO _ '^EM 



can be written as a 
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with K{s) a known function and 
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• as a consequence of CVC, the 7=1 parts 
of both 7?(s) and the hadronic r decay dis- 
tribution provide experimental determina- 
tions of the 7=1 vector spectral function, 
p^^^(s), allowing r decay data to, in princi- 
ple, be incorporated into the determination 

r r -ihad.LO m 

of [a^] [T] 

• after known isospin-breaking (IB) correc- 
tions to the CVC relation have been per- 
formed, the current t and electroproduction 
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(EM) versions of p^~^{s) do not agree [2], 
leading to incompatible determinations of 



The T-based determination corresponds to a SM 
value for compatible, within errors, with the 
BNL E821 experimental result [4], while the most 
recent assessment of the electroproduction (EM) 
determination leads to a value in 3cr disagree- 
ment with experiment 

The discrepancy between the measured EM 
version of p^^^(s) and that implied by IB- 
corrected T decay data is now generated essen- 
tially entirely by contributions from the tt+tt^ 
and TT+TT^Tr'^Tr'' states. The differences, more- 
over, are such that, where the two p^^^{s) de- 
terminations disagree, that based on the r data 
lies higher. Figures of the preliminary BELLE 
T~ — > 7r~7r''z/^ data [51 indicate that, while the 
region over which a discrepancy is observed to 
occur may be shifted by the BELLE results, it 
remains true that, where disagreement exists, it 
is the T version which is higher. 

As pointed out in Ref. [6], the uniformity in 
sign of p:^^^(s) — p^j^f (s) makes it rather easy to 
construct finite energy sum rules (FESR's) which 
distinguish between the two versions of p^^^{s). 
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The form of these FESR's is 
/ w{s)p{s)ds — — (f w{s)Il{s)ds (3) 

Jo 27r J\s\=sa 

where w{s) in an analytic weight function, p{s) = 
p^^^{s), and n(s) is the relevent (charged or neu- 
tral current) 1 = 1 vector current correlator. The 
parameter sq is to be chosen large enough that 
the OPE representation of 11 can be used on the 
RHS of Eq. ([3]). For weights w{y) which are (i) 
functions of the dimensionless variable y = s/sq 
and (ii) both non-negative and monotonically de- 
creasing on the interval < y < 1, it follows that, 
if the T data is correct, then both the slope with 
respect to sq and the magnitude of the EM spec- 
tral integrals for all sq will be too small, relative 
to OPE expectations. Similarly, if the EM data is 
correct, then both the magnitude and slope with 
respect to sq of the r spectral integrals will be too 
large, relative to OPE expectations. It is worth 
stressing that, for the OPE integrals, the slope 
with respect to sq is considerably less sensitive 
than is the normalization to the main OPE input, 
as- We will return to this point in the discussion 
below. 

In the analysis reported here we restrict our 
attention to sq > 2 GeV^, which choice strongly 
suppresses residual OPE breakdown for weights, 
w{s), having a zero at s = sq [7 8\. At these 
scales, the OPE representation is strongly dom- 
inated by the D = perturbative contribution, 
and hence well-determined once is given at 
some particular reference scale. The needed input 
value is obtained by averaging the independent 
high-scale determinations of as{Mz) reported in 
Ref. [9], which yields a^(Afz) = 0.1198 ± 0.0020. 
Values relevant to the lower scales required in 
this analysis are then obtained via 4-loop run- 
ning and matching fTTj. Details of the relevant 
input for, and treatment of, higher dimension 
OPE contributions may be found in Ref. [6] . The 
weights used below, w{y) = 1 — y and u>jv(j/) = 

1 - {wk)y+ (iv^)y'^: N = 3, •••,6, are 
chosen to strongly reduce sensitivity to potential 
poorly known D > 6 OPE contributions [6]. The 
WN^y) also have a double zero at s = sqj fur- 
ther suppressing possible residual OPE breaking 
contributions. 



The 1=1 hadronic t decay distribution has 
been measured by ALEPH pjj, CLEO 13] and 
OPAL [14]. The results below are based on the 
ALEPH determination, for which the covariance 
matrix has been made publicly accessible. Addi- 
tional corrections for long-distance EM effects, as 
evaluated in Refs. |i5j, are applied to the tttt com- 
ponent of the T data [16]. For pre-2003 EM data, 
we employ the results for various exclusive modes 
reported in the compilation of Ref. fT7]. Small 
missing modes are accounted for using isospin re- 
lations and the techniques described in detail for 
such modes in Ref. Post-2003 updates for a 
number of modes |18| have also been taken into 
account. Where uncertainties exist in older publi- 
cations concerning the treatment of radiative cor- 
rections and/or systematic errors, and newer data 
with no such uncertainties are available, we em- 
ploy only the newer data. 

2. Results 

Examples of the results of the comparison be- 
tween the spectral integrals and OPE expecta- 
tions, for the EM case, are shown in Figures [T] 
and [21 for the weights w{y) and wei^y), respec- 
tively. Analogous results for the IB-corrected r 
case are shown in Figures [3| and [4l respectively. 
We see that the normalization and slope of the 
T spectral integrals are in excellent agreement 
with OPE expectations, while both the normal- 
ization and slope are low in the EM case. The 
same behavior is found for FESR's based on other 
non-negative, monotonically decreasing weights, 
though for brevity we have displayed only the 
w{y) and WQ{y) results. As noted above, this pat- 
tern is the one expected if it is the t data which 
are correct and the EM data which are wrong. 
(Note, however, that an alternate r spectral dis- 
tribution which, like the preliminary BELLE r tttt 
distribution, is larger for some s but smaller for 
other s than is the corresponding ALEPH dis- 
tribution, is also capable of satisfying the OPE 
constraints; i.e., the fact that the constraints are 
satisfied by the ALEPH data does not necessarily 
imply that the ALEPH data set is correct.) 

We now quantify the extent of the disagreement 
between the EM spectral integrals and corre- 
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Table 1 

as{Mz) from a fit to the EM and r spectral inte- 
grals at So 4 GeV^ and sq ^ m^, respectively 



w{y) 


EM or T 




w 

W3 

We 


EM 
EM 
EM 


U.iid8_o Q035 

U.ilOZ_Q QQ21 
U.iiOU_Q QQ2g 


1-y 

We 


r 
r 

T 


r, -19-1 9+0. U027 

u.iioy_o.oo2i 

n 11Qc;+0.0020 

U.ii»d_Q QQ22 



Table 2 

The slopes with respect to sq of the EM OPE 
and spectral integrals. The first and second line 
OPE entries for each weight correspond to the 
independent high-scale average and fitted values 
of aa{Mz), respectively, as described in the text. 



w{y) 




SoPE 


w 


.00872 ± .00026 


.00943 ± .00008 






.00934 ± .00008 


We 


.00762 ± .00017 


.00811 ± .00009 






.00805 ± .00009 



Table 3 

Shifts in the effective EM spectral integrals asso- 
ciated with the TTTT and 47r parts of the EM r 
spectral data replacement, for so = 2GeT^^[m^] 



w{y) 


TTTT 


47r 


w 


82% [36%] 


18% [64%] 


We 


87% [45%] 


13% [55%] 



sponding OPE expectations. As a measure of the 
normalization disagreement, we fit as{Mz) to the 
measured EM spectral integrals at so 4 GeV^ 
(the highest scale for which it is still possible to 
construct the EM spectral function as a sum over 
observed exclusive modes) and compare this re- 
sult to the independent high-scale average noted 
above. A similar fit is performed for the r data, 
in this case at the maximum scale, so = m^, ac- 
cessible in hadronic r decay. The results of these 
exercises are shown in Table 1. The fitted values 
are seen to be in good agreement with the high- 
scale average in the r case, but ~ 2 — 2.5 cr low 
in the EM case. Results in agreement with those 
shown in the table are also obtained for other non- 
negative, monotonically decreasing weights ' 19j . 

The results for the EM spectral integral and 
OPE slopes with respect to so {Sexp and Sqpe, 
respectively) are shown in Table 2. In the OPE 
case two values are given for each weight, that in 
the first row corresponding to the independent 
high-scale average as{Mz) = 0.1198 ± 0.0020, 
that in the second row to the fitted value ob- 
tained above for the weight in question, and given 
already in Table 1. We see explicitly, as noted 
above, that the OPE slope is very insensitive to 
as{Mz)- As is evident from the table, no realis- 
tic value for as{Mz) will suffice to bring the OPE 
and spectral integral slopes into agreement; such 
agreement can only be obtained through changes 
in the experimental spectral function. The slope 
discrepancies are at the ~ 2.5 a level. 

To investigate the extent to which the source 
of the normalization and slope disagreement in 
the EM case lies in the / = 1, as opposed to 
7 = 0, portion of the spectral function, we re- 
place the EM TTTT and An data with the corre- 
sponding IB-corrected r results, and rerun the 
OPE/spectral integral comparison. It is found 
that both the slope and normalization of the re- 
sulting r-modificd "EM" spectral integrals are in 
excellent agreement with the OPE constraints. 
An illustration of this point, for the we FESR, 
is given in Figure [5] 

Finally, we consider the relative role of the tttt 
and Att components of the 1=1 r-EM spectral 
function difference in producing the agreement 
between the OPE constraints and the r-modified 



5 



Figure 5. OPE and r-modified "EM" spectral 
integrals for we^y). The solid circles (with error 
bars) represent the original EM spectral integrals, 
the open circles the r-modified "EM" ones. 
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Table 4 

Impact on the EM slope with respect to sq of the 
replacement of EM 47r data with the equivalent 
r data, for w and wq. The modified EM exper- 
imental slope (exp) is given in the first line, the 
OPE slope (OPE) in the second. 



w{y) 


as{Mz) 


Slope 


w 


.1186 


.00936 ± .00026 (exp) 






.00940 ± .00008 (OPE) 


We 


.1176 


.00795 ± .00017 (exp) 






.00808 ± .00009 (OPE) 



"EM" spectral integrals. Since [a^] ' is more 
strongly dominated by the tttt component of the 
spectral function than are the spectral integrals 
appearing on the LHS's of the FESR's employed 
in our analysis, this question is of relevance to 
determining the implications of our results for 

The relative tttt and 47r contributions to the 
spectral integral shifts caused by the pem{s) — + 
Pt{s) replacement are shown, for our two repre- 
sentative weights, in Table 3. The results for the 
other wn are similar to those for wq. The tttt 
component is seen to dominate for sq ~ 2 GeV^ 
and remain important for sq ^ m^. 

Since the agreement between different experi- 
ments for the EM 7r"'"7r~7r°7r° cross-sections is, at 
present, not good [2], it is of relevance to consider 
the impact of replacing only the An component of 
the EM spectral function with the corresponding 
r component. The impact of this change on the 
normalization and slope of the EM spectral in- 
tegrals is shown in Table 4. We see that, were 
future results to bring the EM Att cross-sections 
into agreement with expectation based on r de- 
cay data, the difference between the slope and 
normalization of the EM spectral integrals and 
the OPE constraint values would be reduced to 
the ~ 1 cr level. Such a change would also re- 
duce the discrepancy between the EM-data-based 
SM prediction for and both the T-data-based 
SM prediction and the E821 experimental re- 
sult [20] . The discrepancy between the r and EM 
versions of [a^]'""''^'^ would also be reduced were 
the ~ 4 X 10"^" reduction, relative to the ear- 
lier T average, seen in the preliminary BELLE 
result for the tttt contribution, to remain present 
in the final version of the analysis. Such an im- 
proved EM-T consistency, according to the results 
of this study, would almost certainly be accom- 
panied by a reduction in the difference between 
the experimental value and the EM-based SM 
prediction. Such a development would make even 
more important the role of a reduced experimen- 
tal uncertainty [22 in determining whether or not 
beyond-the-SM contributions have been detected 
in a^,. 
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